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ABSTRACT 

The presence of extremely compact galaxies at z ~ 2 and their subsequent growth in physical size 
has been the cause of much puzzlement. We revisit the question using deep infrared Wide Field 
Camera 3 data to probe the rest-frame optical structure of 935 galaxies selected with 0.4 < z < 2.5 
and stellar masses M, > 10^° '^ Mq in the UKIRT Uhra Deep Survey and GOODS-South fields of 
the CANDELS survey. At each redshift, the most compact sources are those with little or no star 
formation, and the mean size of these systems at fixed stellar mass grows by a factor of 3.5 ± 0.3 over 
this redshift interval. The data are sufficiently deep to identify companions to these hosts whose stellar 
masses are ten times smaller. By searching for these around 404 quiescent hosts within a physical 
annulus 10 kpc < i? < 30 kpc, we estimate the minor merger rate over 0.4 < z < 2. We find 
that 13% — 18% of quiescent hosts have likely physical companions with stellar mass ratios of 0.1 or 
greater. Mergers of these companions will typically increase the host mass by 6% ± 2% per merger 
timescale. We estimate the minimum growth rate necessary to explain the declining abundance of 
compact galaxies. Using a simple model motivated by recent numerical simulations, we then assess 
whether mergers of the faint companions with their hosts are sufficient to explain this minimal rate. 
We find that mergers may explain most of the size evolution observed at z < 1 if a relatively short 
merger timescale is assumed, but the rapid growth seen at higher redshift likely requires additional 
physical processes. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: fundamental parameters — 
galaxies: structure 



1. INTRODUCTION 

The compact nature of massive quiescent galaxies at 
redshifts z ~ 2 was a surprising discovery when it was an- 



nounced some years ago (e.g., 


Daddi et al.|2005l 


Truiillo 


et al. 2006 Buitrago et al. 200S 


van Dokkum et a 


.\2m\. 



have effective radii i?e ~ 1 kpc, 3 — 5 times smaller than 
comparably massive early-type galaxies in the local uni- 
verse. This suggests that they grew significantly in size, 
but much less in stellar mass. Initially there was some 
suspicion that the stellar masses of the z ~ 2 sources 



were overes timated , but deep spectroscopic data ( Cap- 
pellari et al.u2009t [Newman et al.|j2010[ van de Sande 
et al.|20lT fnave verified dynamically the nigh masses of 
selected 1 < z < 2 sources and, in conjunction with the 
abundance of dynamical masses for lower red shift sources 

provided a 



Treu et al. 


2005 


van der Wei et al. 


2005 



valuable, independent confirmation of the size evolution 
Only two physical explanations have been put forward 
to explain this remarkable growth in size while avoid- 
ing the overproduction of present-day high-mass galax- 
ies. Adiabatic expansi on throug h significant mass loss 
can lead to size growth (Fan et al. 2008, 2010). A galaxy 
that loses mass as a result ot winds driven by an active 
nucleus or supernovae, for example, will adjust its size in 
response to the shallower central potential. However, the 
"puffing up" arising from baryonic mass loss occurs only 
when the system is highly active and young in terms of 

|anewman@astro.caltech.edu| 



its stell ar population dRagone-Figueroa &: Granato|2011_ 
see also Bezanson et ai.|2009 ), so it is difficult to see how 
this mechanism can account for the gradual and persis- 
tent growth in size observed for compact sources that are 
mostly quiescent in nature. 

In a hierarchical picture of galaxy formation, merg- 
ers are expected to lead to growth in size and stellar 
mass. Whereas major mergers, involving nearly equal- 
mass components, will lead to comparable growth in both 
size and mass, minor mergers involving lower-mass com- 
panions can produce more efficient size growth (Bezan^ 
20091 |Naab et al.||2009[ [Hopkins et al.||2010c] r 
lanism requires a high rate of occurrence of 



son et al 
I'his mec 

minor mergers, a significant fraction of which must in- 
volve gas-poor companions. Although the major merger 
rate is observati onally constrained reasonably well over 
< z < 1 (e.g.,|Kartaltepe et al.ir2007[ |Lin et al.||2008" 
Bundy et al.|2009[ |de Ravel et al.|2i 



Lotz et al.l 



and via a few measurem ents up to z 3 (e.g., Bluck et al.^ 
2009||Man et al.||2011[ ), the rate at which minor merging 
occurs requires exquisitely deep photometric data. For 
this hypothesis, the key question is whether observations 
confirm that minor merging occurs at the required rate. 

The infrared Wide Field Camera 3 (WFC3/IR) on 
board the Hubble Space Telescope (HST) enables us to 
address the question of whether minor merging is suffi- 
ciently frequent to account for the size growth of com- 
pact sources since z ~ 2. The CANDELS survey (GO 
12444/5; Pis: H. C. Ferguson and S. M. Faber) provides 
an excellent resource for addressing this question since. 



2 



Newman et al. 



in the first two fields to be observed - the UKIRT Ultra 
Deep Survey and southern GOODS fields - the associ- 
ated ground- and space-based and photometry spanning 
0.4 — 8/im is sufficiently deep not only to identify possible 
companions ten times less massive than their hosts, but 
also to reliably determine their photometric redshifts so 
that a physical association can be evaluated. 

Our goal in this paper is thus twofold. First, exploiting 
the unique combination of depth and angular resolution 
in the CANDELS near-infrared data, we aim to measure 
the size growth of massive galaxies. We will show that 
the most compact sources virtually always have quiescent 
stellar populations. We then estimate the minor merger 
fraction by searching for low-mass companions around 
these quiescent sources within a fixed search annulus of 
10 kpc < R < 30 h^^ kpc. A physical association 
can be made through their photometric redshifts. We 
will then interpret the minor merger fraction as a possible 
cause for the growth rate of compact massive galaxies. 

The plan of this paper is as follows. In Section 2, we 
introduce the CANDELS WFC3/IR images and the as- 
sociated photometric data. We describe the selection of 
935 galaxies with stellar masses > lO^'^ '' Mq in the pho- 
tometric redshift range 0.4 < z < 2.5. Section 3 analyzes 
the size growth for this sample and compares our results 
to earlier work. Section 4 introduces our search for faint 
companions around 404 quiescent galaxies spanning the 
redshift range 0.4 < z < 2 in which we can confidently 
detect companions with 10% of the stellar mass of their 
hosts. We discuss the robustness of our search, make cor- 
rections for spurious unassociated pairs, and assess the 
stellar mass content and colors of these companions. Fi- 
nally, in Section 5 we interpret our minor merger rate 
in the context of size growth. After discussing the size 
growth of the quiescent population, we turn to a test 
that asks whether the merger rate is consistent with the 
increasing rarity of compact examples at later times. Fi- 
nally, we summarize our conclusions and the remaining 
uncertainties in Section 6. 

Throughout the paper, we adopt a concordance cos- 
mology with {ilm, ^v,h) = (0.3, 0.7, 0.7 ) and use the AB 
magnitude system ( Oke &: Gunn||1983 ). 



2. DATA AND CATALOGS 

We have compiled an extensive database of optical and 
infrared observations from spac e and the ground in the 
UKIRT Deep S urvey (UD S, Lawrence et al.| [2007l) and 
GOODS-South ( [Giavahsco et al.|2004| ) fields, ottering the 
wide spectral coverage from 0.4 to 8/im necessary to se- 
cure quality photometric redshifts, stellar masses, and 
stellar population parameters for mass-complete samples 
of galaxies to z ~ 2.5. Although our supplementary pho- 
tometry covers a much wider area, we restrict our atten- 
tion to the CANDELS WFC3/IR footprints, since our 
program requires the depth and angular resolution in the 
rest-frame optical afforded by HST. 

2.1. Imaging Data 

The UDS and GOODS-S fields have been observed 
with HST /WFC3 in the J (F125W) and H (F16 0W) fil- 
ters ( |Grogin et al.||2011[ [Koekemoer et al.||2011[) . In the 
UDS, the vO.5 mosaics of the two epochs ot WFC3/IR 
imaging were coadded. For the Advanced Camera for 
Surveys (ACS) F606W and F814W imaging in the UDS, 



we used only the second epoch of observation, since the 
first epoch contained some reduction artifacts at this 
time of this work. The HST imaging was supplemented 
by deep Subaru BVRiz imaging from the Subaru / X MM- 
Newton Deep Survey (SXDS; Furusa wa et al. 2008), us- 
ing the mosaics prepared by Cirasuolo et ai. ij zUlOp , and 
by iiT-band imaging from the UKIDSS UDS Data Release 
6 (DR6). Deep Spitzer Infrared Array Camera (IRAC) 
data from the SpUDS survey (PI: J. S. Dunlop) allows 
us to access the rest-frame near-infrared to z ~ 3. We 
cross-referenced our catalogs to the SpUDS MIPS cata- 
log using a positional tolerance of 1". 

In GOODS-S, we use the first three epochs of 
WFC3/IR imaging in the CANDELS Deep area and the 
first epoch of the Wide region. To this we add the 
GOODS BViz ACS imaging, as well as ground-based 
data in U, R and K from V IMOS ( Nonino et al][2009 | 
and ISAAC ( [Retzlaff et aLpOlO) at the Very Large I'ele- 
scope (VLT). I'he two epochs ot ultradeep IRAC imaging 
from the Spitzer GOODS Legacy Science Program (PI: 
M. Dickinson) were co-added to produce a single mosaic. 
We again cross-referenced our catalog to the MIPS cat- 
alog. 

2.2. Catalogs 

For the main photometric catalog, we chose the WFC3 
H band as the detection image, thereby taking advantage 
of the high-resolution HST imaging while maintaining 
a selection that is as complete in stellar mass as pos- 
sible. The H mosaic, distributed on a 60 mas pixel 
scale, was rebinned to a 120 mas scale, and all other 
imaging was registered to this grid. Object detection 
and photometric precision are insignificantly affected by 
this slightly coarser sampling, but the computational ef- 
ficiency is greatly increased. For measurement of struc- 
tural parameters, where the highest possible resolution 
is critical, we created catalogs for each HST mosaic at 
the original scale (60 mas for WFC3/IR and 30 mas for 
ACS) and matched these to the main catalog. 

Each image (ground, HST, and IRAC) was first regis- 
tered to the H-hand mosaic using smooth transforma- 
tions as determined by the IRAF task geomap. The 
images were then drizzled onto the uniform grid, pre- 
cisely conserving flux, using geotran. A composite point 
spread function (PSF) was constructed in each image 
by stacking suitably normalized cutouts of bright, un- 
resolved sources. Matching PSFs is critical for accu- 
rate colors across images of widely varying resolution, 
yet one wishes to avoid unnecessary degradation of the 
high-resolution data as far as possible. We struck the 
following compromise: the ACS and WFC3 J images, 
each of higher resolution than the detection H image, 
were convolved to match the 77-band PSF. Colors were 
then measured in fi xed apertures of 1"5 dia meter by run- 
ning SExtractor (Bertin &: Arnouts|1996 ) in dual image 
mode. For the lower-resolution imaging from ground- 
based instruments and IRAC, we measured the X band 
flux fx in a wider aperture (see below) appropriate to 
the PSF in a given band X. We then convolved the H 
image to match the X PSF and measured the H flux 
/if.widc in the wide aperture. Finally, fx was scaled by 
the ratio f^ j"^/ fH.widc, in order to refer all fluxes to a 
common aperture. In this way, the HST resolution is 



Can minor merging account for the size growth of quiescent galaxies? 



3 



degraded as minimally as necessary for each band. 

To determine a convolution kernel that matches two 
PSFs, we took the analytic Moffat kernel that best 
matched the curves of growth, weighting toward the radii 
relevant for our aperture photometry. This method typ- 
ically matched curves of growth to ~ 1 — 2%. Colors 
between HST filters were measured in fixed apertures of 
I'.'S diameter. For broader PSFs, the aperture diameter 
was set proportionally to the size of the PSF: 4x the 
half-light radius, but restricted to lie within the range 
I'.'S — 3". The upper limit was chosen to avoid exces- 
sive confusion in the IRAC data. Aperture colors were 
scaled to total fluxes using the SExtractor AUTO aper- 
ture in the iJ-band image. Photometric uncertainties 
were determined using apertures placed at random in 
blank sections of the images. Systematic uncertainties of 
4% (10% in the IRAC bands) were added in quadrature 
to account for zcropoint errors, aperture mismatch, and 
color-dependent flat-field errors in IRAC. Small Galactic 
ex tinction corrections w ere made based on the dust maps 
oflSchlegel et al.|p98|. 



(iTreu et al. 


2010 van Dokkum & Conroy 2010 


Auger 


et al.||2010b 


JNewman et al.||2011 Spinielio et ai. 


;^uii). 



2.3. Photometric redshifts and other derived parameters 
Using this photometry spanning 0.4— 8/^m, pho tomet- 



mer et al. 


2008 


all templa 


;es m 



2008). We permitted linear combinations ot 



based on A'-band flux. Spectroscopic redshift surveys 
have been conducted with the VLT in GOQ DS-S by 



Vanzella et al. ( 2008|^ Popesso et al. ( 2009 ) , and Balestra 



et al.l ( |2UlO| 7while ]Wuyts et al | ( |2b08[ ) have compiled 
redshifts from a number ot additional sources. In the 
UDq^ we draw from C. Simpson et al. (in preparation), 
M. Ak iyama et al. (in preparation), and Small et al. 
(2008). Only spectra with high quality flags were m- 
ciuded. These spectroscopic data provide an opportu- 
nity to test the accuracy of our photometry by forcing 
EAZY to flt templates at the known reds hifts and av- 
eraging the residuals in each filter (e.g., Capak et al. 
[2007). The resulting systematic offsets were small (typ- 
ically < 0.03 mag), confirming the high quality of the 
photometric calibration and PSF matching. The one 
exception was the VIMOS i?-band image , to which we 
added a —0.10 mag correction. In Section [4Tj we assess 
the accuracy of our photometric redshifts by comparing 
to this spectroscopic database. 

Stellar population parameters, including stellar 
masses, were measured by fitting the latest S. Chariot 
& G. Bruzual (2007, private communication) mo dels to 
the broadba nd photometry using the FAST code (Kriek 
et al.|2009b ) . A large grid of models with exponentially 



declining star-formation histories was created, with red- 
shifts between 0.01 and 7 in steps of 0.01(1 -I- z), ages be- 
tween t = 10'' and 10^°'^ years (always less than the age 
of the universe) in 32 logarithmic steps, star-formation 
timescales t between r — 10^ and 10^° years in 31 loga- 
rithmic steps, and dust content varying b etween Ay = 
and 3 in 31 steps. Solar metallicity, the [Calzetti et al.| 
( 2000 ) extinction law, and a Salpeter initial mass func- 
tion (IMF) were adopted. We chose the Salpeter IMF 
because it may be more appropriate for massive galaxies 

^ |http : //www .nott Ingham. ac .uk/astronomy/UDS/data/drS. I 
|html| 



but our analysis is insensitive to this choice since we re 
quire only relative stellar masses. Rather than adopt- 
ing the stellar population parameters of the single best- 
fitting model, we obtain the mean of each parameter by 
marginalizing over the likelihood function. Finally, rest- 
frame colors were computed using the InterRest code 
( [Taylor etal][2009l ). 

I:''igure[T]displays photometry, spectral energy distribu- 
tion (SED) fits, redshift constraints, and color compos- 
ite images for several representative massive galaxies at 
1 < 2 < 2.5. Note that the signal-to-noise ratio is very 
high, even at z ~ 2, reflecting the high quality of the 
photometric data. 

2.4. Survey Mass Limit and Completeness 

We define a limiting stellar mass for our galaxy sample, 
motivated by the desire to obtain a complete census of 
satellites with stellar mass ratios /i* — MgatZ-^^host > 0.1 
at z < 2 as well as our desire to track evolution in the 
sizes of mass-selected hosts to z ~ 2.5. 

The completeness of our catalog was assessed by in- 
serting synthetic objects into blank sections of the UDS 
WFC3 H image, blurring by the empirical PSF and bin- 
ning to the same pixel scale. These were then detected 
using the same SExtractor configuration. The 90% pho- 
tometric completeness limits are iJauto = 26.5 for point 
sources and i?auto — 25.6,25.8, and 26.1 for de Vau- 
couleurs profiles with Re = 0'.'4,0'.'2, and O'.'l, respec- 
tively. For de Vaucouleurs profiles with = O'.'l, which 
is roughly the size expected for local log Af* ~ 9.7 early- 
type galaxies viewed at z ~ 2, the 90% completeness 
limit is i?auto = 26.1. Figures [2^ shows that selecting 
satellites with logAf* > 9.7 at z < 2 ensures iJ-band 
fluxes above this limit, even for a maximally old popu- 
lation. Since we demand completeness for /i* > 0.1, this 
in turn implies a limit of logM* > 10.7 for the hosts. 

If we are only concerned with studies of the host galax- 
ies, i.e., without the need to detect their faint compan- 
ions, they can be followed to somewhat higher redshift. 
We limit ourselves to z < 2.5 in order to retain deep 
detections in F160W, suitable for robust size measure- 
ments at our mass limit. Figure [2]d shows that, in the 
redshift range 2 < z < 2.5, we remain complete at 
log A/* > 10.7 even for Re — 0'.'4, the most extended 
profile we tested. This size corresponds roughly to the 
size of a local = IO^^Mq early-type galaxy viewed at 
z = 2. 

2.5. Surface Photometry and Effective Radii 

We use Galf it ( Peng et al.[[2010 ) to fit Sersic profiles 
to galaxies in our sample, using an automated procedure 
to fit adjacent objects simultaneously. The Sersic index n 
was restricted to 0.5 < n < 8, and the size of the fitting 
box was set by requiring it to enclose the Kron ellipse 
enlarged by a factor of 2.5. The background was mea- 
sured in a rectangular annulus extending 40 pixels from 
the boundary of the fltting box. In order to measure 
structural parameters at similar rest-frame wavelengths, 
we selected different fllters for fltting according to the 
redshift. In the UDS, sizes are measured in F814W for 
0.4 < z < 0.9, F125W for 0.9 < z < 1.8, and F160W 
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Figure 1. Montage of representative massive galaxies at 1 < z < 2.5. High-precision photometry spanning 0.4 — 8fim is plotted along with 
the best-fit spectral synthesis model as described in the text. Composite HST images in the IJH filters (where I is F814W or F775W) are 
inset along with EAZY photometric redshift distributions. 

for 1.8 < z < 2.5. In GOODS-S, F775W is used for 
0.4 < z < 0.75, F850LP for 0.75 < z < 1.1, F125W for 
1.1 < z < 1.8, and F160W for 1.8 < z < 2.5. This en- 
sures that the wavelength at which sizes are measured al- 
ways falls in the rest-frame interval 4240 — 6570 A. Based 
on the mean difference between the Sersic and AUTO 
magnitudes in the H band, we applied slight adjustments 
of AlogA/* = 0.014n to account for light outside of the 
AUTO aperture. 

An extensive suite of tests performed by randomly in- 
serting synthetic Sersic profiles into the F814W, J, and 
iJ-band images showed that we are able to recover radii 
with a t ypical accuracy of 5-10% , consistent with othe r 




90% completeness 

I.I ^r ^ '^r 




studies (van der Wei et al. 



2008 Newman et al. 2010 1 



8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 
log Stellar mass [solar] 

Figure 2. Our sample is designed to ensure a complete census 
of satellites with mass ratios fit, = Msat/Mhost > 0-1 at 2 < 2. 
The relation for a maximally old, dust-free stellar population using 
the Chariot & Bruzual (2007) models is shown as a dashed line, 
while the solid line indicates the 90th percentile in faintness at 
a given stellar mass. These are compared to completeness levels 
(horizontal) to set appropriate stellar mass limits. The top panel 
demonstrates that restricting hosts to log Af» > 10.7 ensures strong 
detections in H for /i» > 0.1 satellites at 2 < 2. The bottom panel 
demonstrates the hosts themselves can be reliably studied to a 
higher redshift of 2 = 2.5. 



This procedure automatically incorporates errors arising 
from background misestimation and blending with neigh- 
boring objects, but applies strictly only to symmetric, 
Sersic-like profiles. In the iJ-band image, we addition- 
ally tested for possible errors caused by PSF variations 
by convolving the synthetic profiles with stellar images 
selected from throughout the mosaic. These were then 
fitted using the empirical stacked PSF used to analyze 
the real data. We found that radii as small as 0'.'05 (0.4 
kpc at z = 2) can be reliably recovered. 

All galaxies with stellar masses exceeding 10^° '' 
were fit. For our study of size evolution presented in 
Section 3, we exclude galaxies for which HST imaging 
in the appropriate filter is not available due to imperfect 
overlap among the observations (5.7% of the sample), as 
well as those whose proximity to the image border or to 
a bright foreground star or galaxy precluded a reliable 
measurement (2.8%). Note that these cuts are uncor- 
related with any galaxy property. We also exclude the 
5.7% of remaining galaxies that are fit with a Sersic index 
n = 0.5 or 8, i.e., the boundaries of the allowed range of 
n. These size measurements are likely to be unreliable. 
Although excluding them may slightly bias our mean size 
measurements, we expect any effect to be minor owing 
to the small fraction of the sample that they represent. 

Effective (half-light) radii are typically reported in a 
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circularized form defined by i?e,circ = a^/?' ^^i^re a is 
the semi-major axis of the half-light ellipse and q = b/a 
is the axis ratio. We adopt a slightly different definition: 
Rh = a{l+q)/2. Physically, Rh closely approximates the 
half-light radius obtained from a classical curve of growth 
analysis on the intrinsic (PSF-deconvolved) Sersic profile, 
i.e., the radius of the circle containing half of the total 
light, as we verified numerically. This definition differs 
appreciably from the more common i?e,circ only for small 
q, for which the latter diverges from a curve of growth 
measurement. For our mass-selected sample, the mean 
(median) difference between Rh and i?e,circ is only 5% 
(2%) and has no impact on the evolutionary trends that 
are the main subject of this paper. 

2.6. Comparison to the Sloan Digital Sky Survey 

The total area covered by our UDS and GOODS-S 
catalogs is 311 arcmin^. At z < 0.4, too little vol- 
ume is probed to provide reasonably large and repre- 
sentative samples of galaxies. In the following analy- 
sis, we therefore supplement our catalogs by compar- 
ing to z ~ galaxies in the Sloan Digital Sky Survey 
(SDSS DR7; [Aljazajian et al.|[2009l ). We selected galax- 
ies from the spectroscopic survey m the redshift interval 
0.05 < z < 0.07. These were matched to stellar mass and 
star for mation rate estim ates from the MPA-JHU DR7 
catalog (Kauff mann et al.||2003[ p] and to Sersic fits from 
the NYU Value Added Catalog ( [Blanton et~al] |2005[ ) . 
The stellar masses were shifted by -1-0.19 dex to convert 
from a Kroupa to a Salpeter IMF. 

There may be substantial systematic differences be- 
tween the derived measurements in the SDSS and CAN- 
DELS. For example, our SED fits include NIR photom- 
etry, while the SDSS do es not. Comparis ons of effective 
radii are also uncertain, 
files to SDSS images o: 
ies. Around lO^M© 



Guo et al. (20091 fit Sersic pro 



representative massive galax- 
their effective radii are on aver- 



age 0.2 dex larger than the Blanton et al. (2005) values. 
Since none of the results in this paper rely on the SDSS 
data, we simpl y adop t the MPA-JHU stellar masses and 
Blanton et al. (2005) radii and, where appropriate, we 
caution how uncertainties in these affect the analysis. 

3. SIZE EVOLUTION OF MASSIVE GALAXIES 

The unique depth, resolution, and area of the CAN- 
DELS near-infrared images provides an opportunity to 
freshly examine the rate of size growth for various cat- 
egories of galaxies within our mass-selected sample over 
0.4 < 2 < 2.5. Below we will focus on evolution in the 
stellar mass - size plane: 



M, 
IQiiM, 







(1) 



In the nomenclature of early-type galaxies, this is the Ko- 
rmendy projection of the stellar mass fund amental plane 
(relating M*, Rf^ and cr; e.g.. Auger et al. 2010a). It 
has been extensively studied, particularly at high red- 



jection (e.g 


., Trujillo et al. 2 


006; Toft et al. 2007; Trujillo 
Buitrago et al.,2008i|CimattT 


et al. 


2007| 


Zirm et al.|ii007, 


et al. 


12008 


Ivan Dokkum et al.||2008| |van der Wei et al. 



20081 |Dam 


anov et al.| 12009 


Toft et al.||2009| |Mancini 


et al. 112010 


Ryan et al.||2010 


Saracco et al. 2011 Dam- 


janov et a 


|2011|). The mass-size plane provides some 



of galaxies (e.g., Nipoti et al.||2003 ), which we exploit in 
Section 5. 

Our sample contains 935 galaxies in the interval 0.4 < 
z < 2.5 with stellar masses exceeding logM* = 10.7. 
Figure [3^ demonstrates a strong correlation between size 
and the specific star-formation rate (SSFR, the star for- 
mation rate per unit stellar mass), such that the most 
compact galaxies are the most quiescent. The lower en- 
velope of points delineates an evolving "compactness" 
limit. This figur e confirms the results of many previ- 
ous studies (e.g.,|Trujillo et al.||2006"Franx et al.||2008 



Williams et al. mm IWeinzirl et al. 2011 but represents 



an important advance, since it is based on a large, ho 
mogeneous sample with space-based sizes uniformly mea- 
sured in the rest-frame optical to z = 2.5. The advantage 
of space-based imaging is particularly evident for lower- 
mass galaxies with logilf* < 11. Most of these that are 
quiescent at z > 1.4 have radii comparable to or smaller 
than (X'l — 0'.'2, which is generally taken as the l imit for 
reliable size measurements in se eing-limited data ( Bezan- 
son et al.||2011| [Williams et al.|[2011| ) 



We expect the appearance of the mass-size plane to 
change with time both through the evolution of existing 
galax ies and the continued emergence of new s ystems 



(e.g., Robertson et al.||2006 Hopkins et ari|2010b). Nev 



ertheless, the evolution ot ihe compactness threshold is 
strong enough that by z ~ 2.5, the most compact galaxies 
are typically smaller than any galaxy found in the lowest 
redshift bin. Although there may be a few compact sys - 
tems persisting even to z = ( Valentinuzzi et aL||2010 ), 
their comoving numb er density is clearly g reatly depleted 



( [Trujillo et al 
that individua 



2009 Taylor et al. 2010). This implies 



compact high-z systems must grow in 



http : //www.mpa-garching.mpg. de/SDSS/DR7/ 



size, and that the responsible processes must evacuate 
the most compact regions of the mass-size plane at a 
rate consistent with Figure [3] For this reason, in the 
following we concentrate foremost on quiescent galaxies, 
which are the most compact. 

Figure [3|d shows the trends we find for 483 quies- 
cent galaxies, defined as the subsample with SSFR < 
0.02 Gyr~^ and no detection in the MIPS 24/xm channel, 
which would indicate the presence of warm dust. Several 
other definitions of quiescence are common in the liter- 
ature. Among these, we note that 88% of our quiescent 
sample w ould be selected by th e UVJ color cuts intro- 
duced by Williams et al. (2010). The median Sersic in- 
dex of the quiescent subsample evolves modestly, from 
(n) ~ 3 to 4.5 over our entire redshift baseline, while the 
median axis ratio is essentially constant at (g) = 0.66. 
This is consistent with the majority of these galaxies be- 
ing bulge-dominated, although some are sur ely disks (see 
Kriek et al.[|2009a[ [vander Wei et ar][20TI| . 

Solid lines in l^'igure [3[ show fits to Equation [TJ which 
are reported in Table[l] Interestingly, there appears to be 
little or no evolution in the slope /3 of the mass-radius re- 
lation within the present uncertainties: formally, we find 
d^/dz = 0.05±0.10. Further, the mean = 0.61±0.05 
is consistent with the /3 = 0.57 we measure for galaxies 
selected in the SDSS using the same stellar mass and 
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Table 1 

Fits of the Mass-Size Relation of Quiescent Galaxies to 
logi?h=7 + /3(IogAf. -11) 



Redshift 


7 


/3 




SDSS z = 0.06 


0.54 


0.57 


0.16 


0.4 < z < 1.0 


0.46 ±0.02 


0.59 ±0.07 


0.21 ± 0.01 


l.0< z < 1.5 


0.30 ± 0.02 


0.62 ± 0.09 


0.23 ±0.02 


1.5 <z < 2.0 


0.21 ± 0.02 


0.63 ±0.11 


0.24 ± 0.02 


2.0 < z < 2.5 


0.04 ± 0.04 


0.69 ±0.17 


0.26 ±0.03 



Figure 3. Size evolution of massive galaxies over 0.4 < z < 2.5. (a) All galaxies with logA/» > 10.7, with color encoding the SSFR. At 
each redshift there is a strong relationship between SSFR and size, with the most quiescent galaxies being the most compact, (b) The 

quiescent subsample, with color now encoding the extinction Ay. Linear fits show the best fit to Rf^ oc with /3 as a free parameter 

i solid line) or fixed to the slope /? = 0.57 (dashed). Dotted lines indicate the la vertical scatter. Spectroscopic samples from [Newman et al.| 
2010[ l and |van Dokkum et al.| | |2008[ l (using CB07 fits from [Muzzin et al.|2009[ ) that pass our selection criteria are plotted as diamonds 
and squares. Sizes represent Sersic effective radii measured at rest-frame ^ 5000 A as described in Section 2. 

SSFR criteria. In the context of spheroids, it is known 
that this slope cannot be es tabhshed solely by dry merg- 
ers of smaller systems (e.g., Ciotti et al.||2007" ), and that 
it must therefore be imprinted by dissipational processes 
during a spheroid's formation, i.e., before it becomes qui- 
escent. From this perspective, it is perhaps expected that 
the mass-radius slope for quiescent systems should per- 
sist to very early epochs. 

Fits to the mass-size relation are always subject to 
an Eddington bias arising from the steep mass function. 
This steepness implies that near the limiting mass thresh- 
old, lower-mass galaxies are scattered above the thresh- 
old more frequently than higher-mass galaxies are scat- 
tered below it. We estimated this bias through Monte 
Carlo simulations, generating mock data with errors in 
stellar masses and radii typical of our sample. These 
were fit to a linear relation using a simple least-squares 
regression with equal weighting, as was done for the real 
data. The measured /3 may underestimate the true slope 
by 0.02 — 0.05. Since this correction is small, sensitive to 
the true errors in the stellar mass estimates, and similar 
at each redshift, we decided not to apply it. 

Noting the lack of significant evolution in the slope 
of the mass-size relation of quiescent galaxies, we fix 
/3 — 0.57 (the SDSS slope) and consider the growth of 
the normalization 7 in Figure |4^. This figure displays 
the mean size of quiescent systems normalized to a stel- 
lar mass of 1O^^M0. It is important to recognize that 
the figure concerns the size evolution of the population 
as a whole and not necessarily the growth rate of any 
individual galaxy. Accordingly, we note that the growth 



Note. — Fits are plotted in Figure |3] Errors are determined 
from bootstrap resampling (negligible in the SDSS). The observed 
scatter is measured using the standard deviation. 



rate at fixed mass dlogj/dt accelerates over this inter- 
val, remaining fairly gradual at z < 1 and then noticeably 
increas ing over z ~ 1 — 2.5. We reached the same conclu- 
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^ ( [2010 )■ Figure |4p shows the same 

data plotted against reashitt; there is no apparent change 
in d\og^/dz. We concentrate here on the evolution per 
unit time because it most directly relates to the effects of 
mergers. The blue points in Figure indicate the sizes 
of the star-forming systems in our mass-limited sample. 
Interestingly, the evolution in size is similar to that for 
the quiescent galaxies, so that star-forming galaxies are 
always, on average, a factor of ~ 2 larger than quiescent 
syst ems of the same mass over the entire redshift range 
(see [Law et al.|[20lT| ) . 

Figure ffi) compares our results on quiescent galaxies 
to several recent studies. Overall, there is a fair de- 
gree of convergence given the diverse nature of the sam- 
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Figure 4. Left: Evolution in the mean size of quiescent (red) and star-forming (blue) galaxies, measured at matched rest-wavelength 
and normalized to Af» = IO^^Mq using the slope Re oc M^'^'' ■ Error bars indicate the lu uncertainty in the mean, accounting for random 
sampling errors only. The shaded region shows the Icr scatter in the qui escent populatio n as measured in Table 1. The large red dot 
indicates our default SDSS relation; the arrow estimates the change if the [Guo et al.||2009| sizes were used instead (Section 2.6). Right: 
Our results (red circles) are compared to other recent estimates, as indicated by the legend. 

pies, which apply various selection techniques to different 
types of data (e.g., sizes measured in different wavebands, 
from space and the ground, selection by color or mor- 
phology). In compiling these data we have harmonized 
all stellar masses to a Salpeter IMF and have applied 
an additional correction of A log M i, — — 0.05z for data 
fit with [Bruzual h Chariot] ( |2003| BC03) models]^ We 
caution that direct comparisons ot simple parametric fits 
may be misleading, since these can depend strongly on 
the redshift interval that is fit. 

The primary conclusion from the high-quality CAN- 
DELS data now in hand is a factor of 3.5 ± 0.3 growth 
in size at fixed stellar mass for quiescent sources over the 
redshift interval 0.4 < z < 2.5, with evidence for acceler- 
ated growth at earlier times (Figure|4^). Our challenge in 
the remainder of the paper will be to attempt to explain 
this growth rate. Although most workers have focused on 
the growth of the mean size at a given epoch (Figures 4), 
there is valuable information in the distribution of sizes 
which can be used to discriminate between the growth 
of individual systems over time and the arrival of new 
members of the population. Although we will discuss 
this model in more detail in Section 5, it is helpful to de- 
scribe the data in terms of the evolving size distribution 
at this juncture. 

Figure Js] shows the cumulative and differential (inset) 
distributions of the mass-normalized radius 7 for quies- 
cent galaxies in several redshift bins. The distribution is 
positively skewed in the higher redshift bins, i.e., it ex- 
hibits an excess of galaxies with large 7, which is mostly 
clearly visible in the inset. The largest quiescent galax- 
ies at a given mass frequently show signs of dust (see 
coloring in Figure [sja), suggesting that their rest-optical 
sizes are impacted by central attenuation]^ For our study 



-0.2 0.0 0.2 0.4 0.6 0.8 
IVIass-normalized radius log y = log Rf, / U° f 



Figure 5. Cumulative distribution of the mass-normalized radius 
7 = Rh/M^^"^ in several redshift bins. The model described in the 
text is overlaid in bands with widths indicating the 90% confidence 
interval. Differential distributions are inset. 

in Section 5, the driving quantity is the declining abun- 
dance of compact galaxies. Therefore, when fitting the 
size distributions, it is important to adopt an asymmet- 
ric form so that the distribution at small 7 is not affected 
by a few apparently large galaxies. 

With this in mind, we describe the size distribution at 
a given redshift with a model in which log 7 follows a 
skew normal distribution. The skew normal distribution 
has three parameters: the mean (log 7), the standard 
deviation criog^, and a shape parameter s that is related 
to the skewness. Appendix A summarizes the relevant 
mathematical details. We parameterize the evolution in 
each parameter as linear in redshift: 



(log 7) (2;) = (log 7) 



d(fog7) 
dz 



(^-1), 



(2) 



^ This accounts for the average difference between BC03 and 
CB07 stellar mass estimates in our quiescent sample. The red- 
shift dependence is expected, since the TP-AGB phase that distin- 
guishes these models is predominant at ages of ^ 1 Gyr. 

* To illustrate the effect of extinction, if we restrict to the ~ 80% 
of quiescent galaxies with Av < 0.6, the intercepts in Table 1 
decrease by A log 7 = —0.05 at ^ > 1, the slopes vary by < Icr, 
and the scatter becomes cri^g^ = 0.20 dex in every redshift bin. 



and similarly for criog-y and s. We then used a Markov 
Chain Monte Carlo procedure to sample the likelihood 
function. Each galaxy was weighted inversely to the 
number of galaxies at similar redshift to ensure that the 
entire redshift range contributed equally to the fit. Fig- 
ure [5] compares the observed distribution in log 7 to the 
model with parameters listed in Table [2] 
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Table 2 

Size Evolution Model 



Mean: (log 7) (2: = 1) 

Standard deviation: criog-y(2:- 

Shape: s(2 = 1) 
dsjdz 



1) 



0.38 ±0.01 
-0.26 ±0.02 
0.22 ±0.01 
0.044 ±0.017 
2.3 ±0.4 
1.0± 1.1 



Note. — Mean quantities, marginalized over all other parame- 
ters, are reported along virith their \a uncertainty. 

This simple model accurately captures the observed 
features of the size evolution. First, (log 7) evolves nearl; 
linearly in redshift as — 0.26z, which Figures and 
demonstrate is a good fit. Second, the scatter criog^ 
evolves fairly little with redshift. The mild increase is 
driven mostly by the increasing abundance of large, dusty 
systems toward higher redshifts, as discussed previously 
(see footnote 4). Note that we have not attempted to 
deconvolve errors arising from uncertainties in the stel- 
lar masses and radii of individual galaxies. Assuming the 
formal stellar mass uncertainties and a 10% uncertainty 
in the radii, the error in individual log 7 measurements 
would be 0.07 dex nearly independent of redshift. Since 
this is much smaller than the measured width of the dis- 
tribution, the intrinsic widths would be only ~ 0.01 dex 
smaller than the measured ones. If the true errors were 
instead twice these estimates, the intrinsic widths would 
be ~ 0.05 dex smaller than the measured ones. The im- 
pact of measurement errors is developed further in Ap- 
pendix B. 



4. SATELLITES OF QUIESCENT GALAXIES AT 0.4 < Z < 2 

The most frequently invoked and well-motivated phys- 
ical process behind the strong, regula r size evolution pre 



sented in Section 3 is merging (e.g., van dcr Wcl et al. 
^0081 IBezanson et al.||2009i [INaab et al.l|2009; Hopkins 
et al.| 2010b). Most previous studies of merger rates 
have been confined to z < 1.4 or have focused on "ma- 
jor" mergers with stellar mass ratios /.t* > 0.25. This 
is partly due to observational limitations, since prob- 
ing higher redshifts and lower-mass companions requires 
deep near-infrared data, and also because major mergers 
are of special interest for studies tracking morphological 
transformations. 



Size grow th, as well as spheroid formation (Bundy 



et al. 



2007 1 , is unlikely to be explai ned by major merg- 



mg alone. Major mergers are rare: Bundy et al. (20091 
estimate a rate of only 0.03-0.08 Gyr"^ for > lO^^-^Mo 
galaxies over 0.4 < z < 1.4. If such low rates persist to 
z = 2, then < 15% of galaxies present at z = 2 will expe- 
rience any major mergers by z = 1, whereas substantial 
size growth must occur over the same period. "Minor" 
mergers involving lower mass ratios may be crucial. 

In this section, we measure the incidence of close com- 
panions to the same set of massive, quiescent galaxies at 
0.4 < z < 2 whose rate of growth was charted in Sec- 
tion 3. As discussed in Section 2, we here limit ourselves 
to z < 2 in order to maintain completeness for stellar 
mass ratios /i, > 0.1. Below, we refer to this quiescent 
sample as the 'primary sample, while the population of 
potential satellites is called the secondary sample. We 
search for secondaries around each primary galaxy at 
projected separations oi 10 < R < 30 h"^ proper 




■(P)' Primary:' z„=1 .7at0.13 
fS) Secondar/: z„=1 .57±0.20 
(I) Interloper: Zp=0.76±0.06 




Wavelength [|jm] 

Figure 6. Demonstration of the pair counting procedure. The left 
panel displays a composite F160W/F125W/F814W image around 
a 1O-'^'''^M0 quiescent "primary" galaxy at Zp = 1.73. The 
10 kpc < _R < 30 kpc search annulus is outlined. One 
/i* ~ 1 : 8 secondary "S" is identified as a possible physical asso- 
ciation based on its consistent photometric redshift (right panel). 
A blue galaxy "I" within the search aperture is excluded based 
on its low photometric redshift. The right panel shows the SEDs 
and best-fitting FAST models. For clarity, the models have been 
smoothed and the fluxes of the interloper reduced by a factor of 
2.5. 

kpc with stellar mass ratios 0.1 < /i* < 1. Note that 
the upper limit avoids double counting. In principle, the 
size of the search annulus should not matter for mea- 
suring merger rates, since the merger timescales increase 
with the search area. In practice, the inner radius avoids 
searching for secondaries buried within the light of the 
primary at low redshift, while the outer radius strikes 
a reasonable balance of finding useful numbers of pairs 
without being dominated by chance alignments. 

Since many galaxies that are close in projection lie at 
different redshifts, we attempt to secure physical associ- 
ations by additionally requiring that secondaries have a 
photometric redshift consistent with the primary, as de- 
tailed below. An example of this is given in Figure [6j 
However, due to the coarseness of photometric redshift 
estimates, some galaxies selected by this method will still 
be chance alignments not physically associated with the 
primary. This contamination rate is estimated simply 
by randomizing the positions of the primaries through- 
out the imaging area, maintaining all their other proper- 
ties, and repeating the search for secondaries using the 
same criteria. This procedure is repeated many times 
to improve the statistical accuracy. Below we distin- 
guish projected secondaries, which comprise all secon- 
daries found within the search apertures, from the sta- 
tistical secondary population that remains after chance 
contaminants are correct for as just described, which we 
term physical secondaries. As we discuss in Section 5, it 
is important to realize that some fraction of these phys- 
ical secondaries will not be bound to their primary host 
and therefore only represent candidate satellites or future 
mergers. 

Below we measure the mean number of physical sec- 
ondaries per primary host and assess the stellar mass 
content and colors of these systems. To examine the 
redshift dependence of these quantities, we break the 
primary sample into three redshift bins of z — 0.4 — 1, 
z = 1 — 1.5, and z = 1.5 — 2. In Section 5, we turn to the 
question of whether the size growth measured in Section 
3 is consistent with the merger rates inferred here. 

4.1. Photometric Redshift Accuracy 
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Figure 7. Photometric redsliift errors in two redsiiift bins. Top: 
Comparison with spectroscopic redshifts for logM* > 9.7 galax- 
ies at 0.4 < 2 < 1 indicating a small scatter cr^z/^i^z) — 
0.024. Bottom: Excess probability of a given redshift difference 
A2:piiot/{l + 2i) for secondaries within 30 kpc of the primary 
quiescent galaxy sample at 1 < z < 2. Here secondaries are se- 
lected based on H flux ratios as described in the text. The quoted 
uncertainty refers to the difference in two photometric redshifts. 
In both panels, the dotted lines indicate the adopted threshold. 



The secondary galaxy sample is selected to have stel- 
lar mass ratios 0.1 < /z* = M2/M1 < 1 and photometric 
redshift differences Sz = {z2 — zi)/{\ + zi) less than a 
fixed threshold, where the subscripts 1 and 2 refer to the 
primary and secondary galaxies. Determining an appro- 
priate threshold for Sz requires knowledge of the accuracy 
of the photometric redshifts. At 0.4 < z < 1, 327 of 1244 
galaxies with logM* > 9.7 (the lowest secondary mass 
we might consider) have spectroscopic redshifts from the 
sources described in Section 2. Figure [7^ compares these 
to photometric redshifts, demonstrating a small scatter 
of crsz — 0.024P] We verified that the redshifts, colors, 
and masses of the spectroscopic subsample at 2; < 1 are 
reasonably representative of the parent population, so 
the measured scatter should reflect the true photometric 
redshift uncertainties. Based on this result, we adopt a 
threshold of \6z\ < 0.1. In 2.8% of cases, the photometric 
estimates differ "catastrophically" by \6z\ > 0.1. 

At 2 > 1 the availability of spectroscopic redshifts 
(^spoc) declines rapidly. Of the massive, quiescent galaxy 
sample, 40 galaxies at z — 1 — 2.3 have measured Zspcc, 
of which only 4 are z > 1.4. The corresponding photo- 
metric redshifts display a small scatter agz = 0.023 with 
only one outlier|^ For the full mass-limited sample with 

^ Throughout, we measure this scatter using the norni alized me- 
dian absolute deviation; see, e.g., |Brammer et al.| l [2b08[ |. 

^ This single 2;spec also disagrees with the photometric redshifts 



logM* > 9.7 at z — 1 — 1.5, the 152 available Zspcc indi- 
cate a scatter of asz = 0.021, while at z = 1.5 — 2 the 32 
available Zgpoc indicate asz = 0.058. We have excluded 
X-ray sources in these comparisons, since they are over- 
represented in the spectroscopic data. We also note that 
the vast majority of z^p^c at 2: > 1 are in GOODS-S, so 
we must assume that similar techniques produce similar 
results in the UDS. Since asz appears to increase toward 
z = 2, we adopt a wider selection \6z\ < 0.2 for select- 
ing secondaries at z = 1 — 2. With this selection, the 
catastrophic error rate {\5z\ > 0.2) is 3 ± 1% and 6 ± 4% 
at z = 1 — 1.5 and 1.5 — 2, respectively, based on the 
available spectroscopic data. 

Since the spectroscopic samples are not representative 
of the full massive galaxy population at z > 1, it is use- 
ful to assess the accuracy of photometric redshifts by 
ot her means. We use the em pirical technique proposed 
by Quadri Sz Williams| (20101. Their method is an ap- 
plication of the general procedure employed throughout 
this section: determine the distribution of 6z for well- 
defined primary and secondary samples, and subtract the 
distribution obtained with scrambled galaxy positions. 
In this situation, it is preferable to define a secondary 
sample based on fiux rather than stellar mass, since er- 
rors in Zp and stellar mass are correlated. To determine 
a limiting flux ratio that best mimics a mass-based se- 
lection 0.1 < /i* < 1, we examined the distribution of 
AH = H2 — Hi between the primary quiescent sample 
and physical secondaries selected based on their stellar 
mass. In 90% of cases, AH < 2.2 mag. This motivates 
a secondary sample defined by < AH < 2.2 mag. Fig- 
ure [TfD shows the distribution of redshift differences 6z 
for the physical secondaries. 

The distribution is broader than at z < 1, as the spec- 
troscopic comparison also indicated. The uncertainty 
f^Sz = 0.047 measured here refers to that in the difference 
between two photometric redshifts. The more impor- 
tant uncertainty for this study is the rate of catastrophic 
{\5z\ > 0.2) redshift errors. A crude estimate of this can 
be obtained by integrating the curve in Figure [Tja, which 
yields 9 ± 15% over z = 1 — 2. Using the same technique 
at z ~ 2, we find a possibly higher catastrophic rate of 
15 ± 20%, but this cannot be determined precisely with 
the present sample size. These noisier estimates may be 
higher than the 3 — 6% inferred from the spectroscopic 
database, but that sample is biased toward bright sys- 
tems. A better assessment of the catastrophic rate will 
require spectroscopic redshifts for larger and more repre- 
sentative samples of galaxies at z = 1—2 than is currently 
available. 

4.2. Subtraction of Host Light 

A concern in all pair studies is that the photometry 
of the secondary galaxies may be contaminated by light 
from the hosts. By inserting synthetic pairs of galaxies 
with 1 : 10 luminosity ratios and projected separations 
10 — 30 h^^ kpc into the iJ-band mosaic, we found that 
our detection efficiency is not affected by the proximity of 
host. Further, these tests indicated that the aperture col- 
ors are less affected than the H-h&nd AUTO magnitude 
used to scale the total stellar masses. To correct for this. 



in the MUS YC jCardamone et al.|2010 [l and FIREWORKS l |Wuyts| 
|et al.|2008| l catalogs^ 
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we measure the SExtractor AUTO magnitudes of the 
secondary galaxies in images from which the light of the 
primary galaxy has been subtracted using our Sersic fits 
(Section 2.5 1. We also compute stellar mass ratios /i* 
using fits that omit the IRAC photometry, which is the 
most susceptible to contamination, although this has lit- 
tle effect on our results. 

4.3. Abundance and Stellar Masses of Physical 
Secondaries 

We now turn to the frequency of physical secondaries 
and their stellar mass content. First, we consider the 
pair fraction /pair- This is simply the mean number 
of physical secondaries per primary galaxy: /pair = 
{Np — Nr)/Ng, where Np is the number of projected 
secondaries, Nr is the expected number of chance align- 
ments given the total search area, and Ng is the number 
of primary galaxies. Throughout our pair analysis, we 
exclude the shallower "Wide" section of GOODS-S and 
primaries for which more than 20% of the search annu- 
lus is masked (e.g., near the image edge). All results 
in the remainder of this section pertain to mass ratios 
0.1 < < 1. 

Table [3| presents the results. For quiescent primaries, 
we find Jpair = 16% ± 3% when averaged over the en- 
tire host mass and redshift range. Moreover, the pair 
fraction does not appear to evolve significantly with red- 
shift within our uncertainties: formally, we find /pair oc 
(1- 



-0.11±0.68 



The paradoxical result that the galaxy 
merger rate remains flat as the halo merger rate in- 
creases with red shift has been explored in many theoret- 
ical works (e.g., Berrier et al.||2006 Kitzbichler & White 
[2008| ). Although our present sample is not large enough 
to be divided in both redshift and mass, we can examine 
possible mass-dependent trends by dividing the sample 
into the three mass bins listed in Table [3] and averaging 
over the full redshift range. We find that the pair fraction 
increases slightly w ith stellar mas s as /pa ir oc jVf^-^^^"'*^, 



in agreement with Bundy et al. (2009). For later use 
in our models of size growth, we also tabulate the "in- 
trasample" fraction /is of physical secondaries which are 
also members of the primary sample (i.e., are quiescent 
and > IQio '^AfQ). 

From the point of view of galaxy assembly, an equally 
useful quantity is the amount of stellar mass contained in 
physical companions. We estimate this simply by com- 
puting the mean total stellar mass in projected secon- 
daries, expressed as a fraction of the host mass, and 
subtracting the random contribution as described pre- 
viously. We denote this quantity Jm. Averaged over all 
masses and redshifts, we find /m — 0.060 ± 0.011. The 
mean mass ratio (/i*) = /j\///pair is very nearly constant 
at « 0.39 in all redshift and primary mass ranges. As 
many authors have noted, this implies that the stellar 
mass delivered in mer gers arrives primarily i n more mas- 
sive secondaries (e.g., Hopkins et al. 2010a). In Section 
5, we compare (yLt*) to theoretical expectations. 

Although we concentrate on the growth of quiescent 
galaxies for the remainder of this paper, for compari- 
son with future work we also tabulate the corresponding 
quantities for star-forming galaxies in Table [3] 

4.4. Colors of Physical Secondaries 



Since physical secondaries likely represent the "build- 
ing blocks" for the future mass assembly of quiescent 
galaxies, particularly in their outer regions, it is inter- 
esting to consider their stellar populations in relation to 
those of their hosts. In particular, the fraction of mergers 
which are "dry" (gas-poor) versus "wet" is an important 
input to models of galaxy evolution. Table [3] presents the 
fraction /q of physical secondaries which are quiescent. 
We calculate this using two definitions of quiescence: the 
SSFR < 0.02 Gyr"^ threshold used throughout this pa- 
per (also excluding MIPS detections), and a color selec- 
tion {U - V)cor > l-l- Here 



{U - y)cor - ([/ - 1^)re.t " 0.47Ay 



(3) 



represents the extinctio n- corrected rest-frame U—V color 
(Brammer et al. 2009). Overall, the two selections are 
qualitatively consistent: most physical companions to 
quiescent galaxies are themselves quiescent at z < 1, and 
this fraction decreases with redshift. 

This is illustrated in Figure |8j which shows the color 
distribution of the physical secondaries. (As throughout, 
we have subtracted the color distribution of similarly- 
selected galaxies in randomly placed apertures.) The 
secondaries are compared to a field sample with matched 
distributions in stellar mass and redshift (solid line) and 
to the primary quiescent host sample (dotted). In both 
redshift bins, the physical secondaries are on average red- 
der than the field comparison sample. The fraction of 
blue secondaries increases with redshift, suggesting that 
the reservoir of future merger candidates includes pro- 
gressively more gas-rich galaxies at earlier times. The 
implications for the merger descendants are interesting 
but not completely clear. On the one hand, secondary 
bursts of star formation are observed in spectroscopi c 
samples of early- type galaxies at z ~ 1 ( Treu et al.|2005 ). 
On the other hand, as we review in section 5, merger 
timescales are expected to be > 1 Gyr. If the processes 
driving satellite quenching are mostly confined to the fi- 
nal ~Gyr, many of these blue secondaries may be much 
redder by the time of the final merger. Nevertheless, it 
seems likely that a significant fraction of mergers at z > 1 
are not completely dry, even for red hosts. 

4.5. Comparison with Previous Work 

Comparisons to independent estimates of the pair frac- 
tion are complicated by the intrinsic differences in sam- 
ples selected by various means (stellar mass, color, lumi- 
nosity). In particular, as we discuss below, samples in 
which satellites are selected based on their stellar mass 
will systematically differ from those based on luminos- 
ity, particularly in the rest-frame optical. An advan- 
tage of the present study is the characterization of the 
merger rate and size growth using a uniform mass-based 
selection. Nevertheless, it is valuable to compare our 
/pair measurements to other works. In the following, we 
rescale published /pair measurements to our s earch area 
by ass uming that /aair(^ < -f^max) oc flrnax ( Kitzbich- 
ler fc White 2008} |L6pez-Sanjuan et al.|[201l| ) i''igure|9 
shows this comparison, tocusmg primarily on those stu(> 
ies that adopted a mass-based selection, included minor 
mergers, or probed to z ~ 2. For comparison, we also 
plot the "major" pair fraction in our sample, defined by 
0.25 < /z* < 1. 
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Table 3 

Abundance and Properties of Physical Secondaries with 0.1 < /x, < 1 



Redshift 


Primary 


Np 


Nr 


^9 


/pair(%) 


/m(%) 


(m.> /q{%) 


/is(%) 




mass range 










= /M//pair SSFR / {U - V)cor 





Quiescent galaxies 



0.4 < 2 < 1 


10.7 < 


log 


M, 


41 


9.7 


177 


18 ±4 


6.2 ± 1.6 


0.35 ± 0.06 


71 ± 10/91 ± 8 


35 ± 12 


1.0 < 2 < 1.5 


10.7 < 


log 


M, 


27 


12.4 


117 


13 ± 5 


5.3 ± 1.7 


0.43 ± 0.08 


66 ± 16/78 ± 18 


33± 17 


1.5 < 2 < 2 


10.7 < 


log 


M, 


30 


11.8 


100 


18 ±6 


7.5 ±2.7 


0.41 ± 0.06 


33 ± 15/38 ± 16 


18± 14 


0.4 < 2 < 2 


10.7 < 


log 


M, < 10.9 


39 


18.0 


141 


15 ± 5 


5.8 ± 1.8 


0.39 ± 0.07 


52 ± 14/71 ± 16 




0.4 < 2 < 2 


10.9 < 


log 


M, < 11.2 


48 


16.5 


174 


18 ± 5 


6.5 ± 1.8 


0.36 ± 0.05 


47 ± 12/67 ± 13 




0.4 < 2 < 2 


11.2 < 


log 


M, 


22 


4.9 


80 


21 ± 6 


8.6 ±2.9 


0.40 ± 0.08 


76 ± 13/82 ± 12 




Star-forming galaxies 




















0.4 < 2 < 1 


10.7 < 


log 


M, 


15 


4.3 


84 


13 ± 5 


4.9 ±2.2 


0.39 ± 0.07 


50 ±23/71 ±21 




1.0 < 2 < 1.5 


10.7 < 


log 


M, 


30 


14.5 


127 


12 ± 5 


4.1 ± 2.0 


0.34 ± 0.08 


41 ±20/53 ±21 




1.5 < 2 < 2 


10.7 < 


log 


M, 


30 


12.9 


108 


16 ±6 


5.7 ± 2.4 


0.36 ± 0.08 


29 ± 15/59 ± 17 





Note. — Np, Nr, and Ng are the number of observed projected pairs, the expected number of these that are chance alignments, and 
the number of primary galaxies, respectively, /pair = {Np — Nr)/Ng is the number of physical secondaries per primary galaxy, of which 
a fraction /q are quiescent (as determined using the two methods described in the text) and a fraction /is are included in the primary 
sample (i.e., quiescent and massive), fuj is the mean stellar mass in physical secondaries as a fraction of the host. Uncertainties in /pair 
reflect Poisson noise in Np and A^,-; for other quantities, uncertainties are determined from bootstrap resampling of the primary galaxies. 
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Figure 8. Distribution of rest-frame colors of physical secondaries around the massive, quiescent galaxy sample in two redshift bins, 
compared to a field sample that is matched in stellar mass and redshift (solid) and to the primary sample (dotted). At higher redshifts, a 
significant fraction of companions are blue. 



At z < 1, our measurements are broadly in agreement 
with previous results when analogous samples a re com- 
pared and search apertures are matched (e.g. , Kartal- 
tepe et al.| f2007^ IL in et al.|[2008l |Rawat et ai:]|2008| |de 
Ravel et al. 2009) . Of particular inte rest is the com 



parison to Lopez-Sanjuan et al. ( 2011 1, who identified 
minor spectroscopic pairs. Figure|^]shows results for 
their red host sample with secondaries having rest-S lu- 
minosity ratios /is > 0.1. A. Nierenberg et al. (2012, in 



preparation, see also Nierenberg et al 



2011) identify mi- 



nor companions to spheroidal hosts split into two bins of 
stellar mass (Chabrier IMF), identifying satellites with 
flux ratios fipsuw > 0.1. The strong mass dependence 
they find highlights the importance of matching hosts in 
stellar mass when comparing pair fractions or analyzing 
the associated size growth. We note that Nierenberg et 
al. use a local background estim ation that is expecte d to 

(|2009| se- 



yield smaller raw pair fractions. Bundy et al. 
lected major pairs at z < 1.4 having K - band nux ratios 
UK > 0.25. Their results for red hosts with log > 10.5 
(Chabrier IMF) are shown in Figure |9] Several authors 
have inferred merger rates from morphological signatures 



(e.g., Lotz et al. 2008a', 'Conselice et al.'2009 Bridge et aLj 
2010|). i^'or a recent review, we refer to, Lotz et al. (2011j). 



"Few other studies have considered minor mergers at 
z > 1. Among these, Williams et al. (2011) is the most 
directly comparable to our work, as their selection is 
based on stellar mass. Figure M shows their quiescent, 
log A/* > 10.8 (Kroupa IMF) host galaxy sample, for 
which they find /pair ~ 0.16 — 0.20 (/i* > 0.1) essentially 
independent of redshift, in encouraging agreement with 
our results. 



Man et al.| pOTT| ) use if-band HS T /NICMOS imag- 
ing to assess the major pair fraction for massive (> 



10 Mq) hosts, selecting se condaries wi th _ff-band flux 
ratios /ig > 0.25. Recently, Bluck et al. ( 2011( ) (see also 



Bluck et al. 2009) studied fainter companions around a 
similarly massive population by identifying close pairs 
to a limiting flux ratio of piH = 0.01 in NICMOS imag- 
ing. As Figure [9] shows, our data are consistent with 
these flux-based selections at z < 1.7. There is a hint 
that the major and total (/i > 0.1) pair fractions rise 
toward z cri 2, but this is not very significant at present. 
( Bluck et al.|2011 demonstrate stronger increases toward 
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Figure 9. Major {fj, > 0.25, small red symbols) and total {fi > 0.1, 
large blue symbols) pair fractions are compared to recent indepen- 
dent measurements at z < 2. Here /i refers to either a flux or 
stellar mass ratio. Overall there is reasonable consistency within 
the statistical uncertainties, with no clear sign of strong evolution 
in /pair over z = 0.5 — 2. Expected systematic differences arising 
from different selection techniques are discussed in the text. Points 
are slightly offset in redshift for clarity. 

z ~ 3.) Furthermore, samples in which secondaries are 
selected based on rest-optical flux wi ll not agree in detai l 
with stellar mass-based samples (see Bundy et al.||2004 ). 
At z ~ 2, the H band probes the rest-trame V band. 
A significant dispersion in the stellar mass-to-light ra- 
tio IvU/L is thus expected, and the mean Af*/L declines 
substantially with decreasing mass. Assuming a constant 
Af* / L equal to that of the host and a limiting flux ratio 
> 0.1, for example, will include bluer galaxies with 
lower mass ratios /x* < 0.1, likely resulting in an elevated 
/pair compared to a mass-selected sample. This effect is 
expected to become stronger toward lower mass ratios 
and toward higher redshifts as the H band probes bluer 
rest wavelengths, possibly impacting trends with redshift 
in flux-sel ected samples 



Finally, van Dokkum et al. (2010) used a novel method 
to infer indirectly the rate at whicn massive galaxies as- 
semble mass through mergers. They tracked the stellar 
mass growth of a sample over z = — 2 with constant 
comoving number density and subtracted an estimate of 
the in situ star formation. Their estimated "specific as- 
sembly rate" is M*/M* = 0.03(1 + z) Gyr"!. At z - 1 
this compares well with our pair counting estimate of 
fu/Te ~ 0.07/re for merger timescales Tg ~ 1 Gyr (see 
Section 5.2), although we find a weaker redshift depen- 
dence. 

5. CONNECTING SIZE GROWTH WITH MERGERS 

In this section we present simple models that compare 
the rate of size growth measured in Section 3 with that 
attributable to mergers of close pairs as studied in Sec- 
tion 4. Sections 5.1 and 5.2 review the theoretical in- 
gredients necessary to convert pair fractions into growth 
rates in mass and size. In Section 5.3 we examine the 
growth in the mean size of quiescent galaxies and dis- 
cuss how the rate of size evolution experienced by any 
individual galaxy may be substantially smaller. Finally, 
in Section 5.4 we combine our constraints on the distri- 
bution of sizes of quiescent galaxies with the evolution 
of their number density to establish a minimum growth 



rate, which we then compare to a merger model. 

5.1. Merger Timescales and the Distribution of Mass 

Ratios 

Converting the observed number of physical compan- 
ions into a merger rate requires us to specify the timescale 
during which a merger appears within our search aper- 
ture, i.e., a projected separation between 10 and 30 
kpc. We deflne an effective timescale that incorpo- 
rates two physical effects. The first is the mean time T^g 
during which a bound, sinking satellite appears within 
our search aperture. As discussed in Section 4, however, 
not all of the physical secondaries we counted are neces- 
sarily bound to their host. By subtracting the number of 
pairs found in randomly placed apertures, we account for 
interlopers in the far foreground and background of the 
galaxies in our primary sample, but we can expect the 
remaining "physical" secondaries to include both bound 
satellites and other galaxies in the larger group-scale en- 
vironment that are not bound. As is common practice 
in merger rate studies, we account for thi s by defining 
a factor Cmg (see also Bundy et al. 20091 to represent 
the fraction of physical secondaries ttiat are bound and 
due to merge on a typical timescale Tmg. The effective 
timescale is then Tr = Tme/ C„ 



.'mg- 



Patton & Atfield (2008) study projected pairs of sim- 
ilar luminosity in the SDSS. They assume that merg- 
ers of luminous pairs occur with a typical timescale of 
Ting = .5 Gyr. Based on test s using the Millennium sim- 
ulation (Springel et al. 20051, they estimate Cme « 0.5 



(their /sd] tor the most iurninous pairs, resulting in an 



effective timescale Te = T^s/Cj 



Lotz et al. (2008b 2010b a 



n |2010b|c 

laxies usn 



1.0 Gyr. 

investigate 



merger 



timescales tor disk galaxies using high-resolution hy 
drodynamical simulations. For self-similar mergers of 
their most massive disk G3, they find a mean timescale 
(Tmg) — 0.7 Gyr within ou r adop t ed search annulus 
Similar to 



(2008), Lotz et al. 



Patton fc Atfield] 
allow for projection effects by setting G'mg 
ing in an effect ive Tg = 1.2 G yr. 



(120111 
resuit- 



Kitzbichlcr fc White] ( |2008[ ) calibrated Te using the Mil- 
lennium simulation, coupled with a semi-a nalytic model 
(SAM) of ga laxy merging and evolution (De Lucia fc 
Blaizot]]2007 |. Considering major mergers (0.25 < /z* < 
4yoF~M* = lO^^iW© galaxies and a search aperture of 
i? < 30 kpc, they find Te = 2.7 Gyr. Since Tg scales 
approximately as the outer radius of the aperture, and we 
exclude the inner 10 kpc, the appropriate timescale 
for our study would be ~ 2/3 of this, or Tg 2 Gyr. 

Taken together, these studies imply Tg = 1 — 2 Gyr 
for major mergers. In order to make progress in the 
present study, we also need an estimate of Tg for minor 
mergers, which will likely be larger. The outcome will 
depend on the relative contribution of minor mergers, 
i.e., the distribution of /i* in our sample, and on how 
strongly the timescales vary with . Figure [lO] breaks 
the observed pair fractions from Section 4 into several 
bins of the mass ratio /i* . The distribution is essentially 
flat in log /u* ; this corresponds to a mass function that 
rises as We note that a uniform distribution over 

— 1 < log^* < has a mean (^*) = 0.39, in agreement 
with our measurement of (/x*) — 0.39 ± 0.04 (Section 4). 
To understand the flat distribution in Figure |10[ we 
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Figure 10. Distribution of the stellar mass ratio ^» of physical 
secondaries around log Af* > 10.7 quiescent galaxies in two redshift 
bins. Merger rate predictions of the^Hopkins ct al. (2010a| SAM are 
overlaid under several assumptions tor tfisTtimescale Tg-. a constant 
Te = 1.2 Gyr (solid line), Te oc fiZ^'^ for fJ-t < 0.3, in line with the 
Lotz simulations (dashed), and Te oc [/.tlog(l + (dotted) 
as expected from pure dynamical friction considerations. 

compare the measured pair fractions t o merger rates 



R 



fpair/Te predicted by the SAM of Hopkins et al 



(2010a I under three choices of the observability timescale 
Te described in the caption. Within the uncertainties, our 
observations are consistent with any of these timescale 
scalings. At higher redshifts, the SAM predicts higher 
/pair (blue curve), which we do not observe. However, 
the distribution in /i^, maintains the same shape nearly 
independ ently of mass o r redshi ft. This is not particu- 
lar to the Hopkins et al. (2()10a [) model, but is a generic 



feature of many SAMs (see Lotz et al.|2011|. The SAM 



also provides an estimate of the amount of stellar mass in 
/z* < 0.1 mergers that we do not probe observationally. 
Only ~ 9% of the predicted mass assembly rate (i.e., the 
rate at which stellar mass is delivered through mergers) 
is due to ^* < 0.1 mergers. This simply reflects the fact 
that for such low mass ratios, the time for a galaxy to de- 
scend from the virial radius to the center quickly exceeds 



a Hub ble time ( [Taffoni et 3172003 ^ 

2008| . By observing > 0.1 pairs, therefore, we expect 



Boylan-Kolchin et al. 



to account for the vast majority of the mass assembly. 

In summary, the effective timescale for major mergers 
is likely Te = 1 — 2 Gyr. For lower mass ratios, estimates 
are even less certain. However, since the physical secon- 
daries are not overly dominated by the lowest mass ratio 
systems, consistent with the predictions of SAMs, we ex- 
pect the appropriate average for our sample to be only 
moderately higher. In the following analysis we present 
results for models spanning a range of timescales. 

5.2. Size Growth Efficiency 

In order to address whether mergers drive the observed 
size growth, we need to know how the half-light radius Rh 
of a galaxy changes after undergoing a merger of mass ra- 
tio ^. This question has been addressed in the literature 
both analytically and using extensive suites of merger 
simulations. The growth efficiency is commonly param- 
eterized by a = d log /d log . For 1 : 1 mergers of 
spheroids, both the mass and rad ius approxi mately dou- 
ble and a ~ 1 (e.g., Hernquist et al.^^1993, ,Nipoti et al 



2003| [Boylan-Kolchin et al [2006^ . 

Simple virial arguments based on energy conservation 
show that the growth efficiency can be higher for more 



minor mergers ( Hopkins et al. 2009 Bezanson et al. 2009 
Naab et al.||200^ ). Assuming that the orbit is parabolic 



that the progenitors and merger product are structurally 
homologous, and that there is negligible energy transfer 
from the stars to the dark halo. 



log(l + ^^-'^) 
log(l + fl) 



(4) 



where we have assumed the progenitors lie on a Rh oc Aff 
relation. For self-similar mergers /i = 1 and we recover 
a = 1. For a mass-radius slope of /3 = 0.57 (Section 
3) and the lowest mass ratios we observe (/i = 0.1), this 
estimate becomes a = 1.6. We therefore expect an ap- 
propriately averaged (a) over the mass ratios we consider 
to lie in this range. 

Recognizing the assumptions entering this simple for- 
mula, it is es sential to verify its p redictions with merger 



simulations. Nipoti et al. (2009) simulated hierarchies 



of multiple dry minor merger s of sp h eroids and found 
(a) = 1.30. More recently, Nipoti (2011) performed 
a suite of = 0.2 dry spheroid mergers and found 
(a) = 1.60 ( see a lso C. Nipoti et al 2011, submitted). 



Oser et al. (2011) investigated the relevance of Equa- 



tion H] in a cosmological hydrodynamical simulation and 
found it to be accurate. Altogether, based on these re- 
sults, we consider a ~ 1.3—1.6 to be a reasonable average 
over the mass ratios we consider. We note that a higher 
efficiency {a > 1.6) has not been demonstrated, when av- 
eraged over a representative set of orbits, in any iV-body 
simulation of which we are aware. 

5.3. Matching the Observed Growth of the Quiescent 
Population 

With estimates for the growth efficiency a and merger 
timescale Tg in hand, we can now proceed to a simple 
model that estimates the rates of growth in mass and 
size due to mergers. In a time interval At, the stellar 
mass of the average quiescent host in our sample in- 
creases by A log M* = log(l -I- /m)^*^^° , while the radius 
by definition increases by a A log Af* . Since we expect 
a > /3 = 0.57, as discussed in Section 5.1, mergers will 
shift the mean mass-size relation: 

A log 7 = A log Rh - /3A log Af, = (a - 0.57) A log 

-(a-0.57)log(l + /M)^*/"= (5) 

We have neglected here the small change in the number 
density arising from mergers within the sample over the 
interval At. This incurs a fractional error of ~ /is/pair ~ 
5% in the mass accreted, negligible for our purposes. 

Figure [TT^ reproduces the observed growth in 7 from 
Figure |4] Using Equation [5j we overlay growth trajecto- 
ries for representative values of the growth efficiency and 
merger timescale to illustrate the evolution of the qui- 
escent galaxy populations in place at z = 2 and z — 1. 
Here we have taken f]\f appropriate to z = 0.4 — 1 and 
z = 1 — 2 (Table l3| and applied a 15% correction to ac- 
count for both aoditional satellites below our /i* = 0.1 
limit an d for possible catastrophic redshift errors (Sec- 
tion 4.1 ). We assume that all galaxies grow smoothly at 
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Figure 11. Left: Observed evolution in size measured at fixed mass for quiescent galaxies, reproduced from Figure |4] compared to 
predicted trajectories expected from the simple merging model discussed in the text. We begin with progenitors at z = 2 and 1 and predict 
their future growth. The width of the shaded bands corresponds to the uncertainty in /jvf . At z < 1, mergers plausibly account for much of 
the size growth if the timescale is short, but at higher redshifts they cannot match the rapid growth of the population. Right: Comoving 
number density of log Af» > 10.7 quiescent galaxies in CANDELS and the NMBS. Error bars for individual fields (slightly offset in redshift 
for clarity) reflect Poisson noise only, while the error in the mean (diamonds) is determined empirically from the scatter among the fields. 
The solid line and shaded region indicate the fit used in our model and its la uncertainty. 



this rate, i.e., we do not incorporate stochasticity in the 

incidence of mergers. 

The primary conclusions from Figure [TT^ are twofold. 
First, at ^ < 1 the pairs we observe can plausibly ac- 
count for most of the observed size growth if an effective 
timescale Tg ~ 1 Gyr, at the short end of the estimates 



There is observational support at z ~ for the idea 
that younger early-type galaxies are larger at fixed mass 



discussed in Section 5.1 and an average growth efficiency 
a w 1.6 are valid. Second, at z > 1 the observed growth 
in log 7 per unit time increases significantly. This en- 
hanced growth rate cannot be matched by mergers using 
any reasonable choices of Te and a. 

As discussed in Section 3, however, an im por tant ob- 
jection to the model comparisons in Figure |11^ is that 
we are tracking the mean growth rate of the entire pop- 
ulation, as if all sources are enlarged in lockstep. In re- 
ality, the population at any redshift comprises both old 
galaxies which formed at higher redshift and which pre- 
sumably are growing via mergers, and sources newly ar- 
riving on the quiescent sequence, whose size may reflect 
their epoch of formation. Galaxies appearing later are 
typically formed from gas-poorer progenitors. They are 
therefore expected theoretically to experience less dissi 
pation in their formation, possibly leading to less com 
pact remnants (Robertson ct al. 2006; Khochfar & Silk 



(jShankar & Bernardi] |2009[ [van der Wei et al.||2009 
Bernardi ct al. 201^7^ Un the other hand, some recent 
studies at higher redshift have found no sign of su ch a 
cor relation (|i>u.Ullo et al. | 2011 [ [Whitaker et al.|2012[ but 



Saraccoet al.||2011[ ). Although the true situation re 



20061 iHopkins et al. 2010b, Shankar ct al. 20lTt: 



l''igure |ll| 3 deinonstrates that the comoving number 
density of quiescent galaxies increases very rapidly at 
z > 1.3, exactly where the growth in mean size is most 
rapid. For example, only ^ 25% of the sample at z ~ 
1 was already formed and quiescent at z 2. These 
early galaxies may need only to grow marginally into 
the compact tail of the distribution at z ~ 1. They 
might then experience significantly less growth than the 
population mean tracked in FigureflTk. In this figure, we 
have combined our CANDELS catalog with t hose from 
the NEWFIRM Medium Band Survey (NMBS; [Whitake7 
[et al.^^2011j ) to increase the total volume. Densities m 
the various fields agree closely at z > 1.5, where large 
volumes are probed, while cosmic variance dominates at 
z < 1.5. 



see^ 

mains unclear, it is interesting to consider size growth 
assuming that the oldest galaxies at a given redshift and 
stellar mass are the smallest, since this corresponds to 
the minimum rate of growth that individual old galaxies 
must undergo. We now seek to construct a test that ac- 
counts for the continual emergence of quiescent systems. 

5.4. A Minimum Rate of Growth for Early Compact 
Galaxies 

The physical processes that determine the size of a 
galaxy in its early history might therefore be quite dif- 
ferent from those that drive its subsequent growth. |Oser| 
et al. (2010) described a "two phase" picture that, while 
obviously a simplification at some level, still provides a 
useful paradigm for galaxy growth. The first phase is 
characterized mainly by in situ star formation, while in 
the second phase, most growth occurs through accretion 
of stars. We wish to test whether mergers are sufficient 
to power size growth in this second phase. As we dis- 
cussed in Section 5.3, growt h in the mean size of the 
quiescent population (Figure 



111) entails processes op- 



erating in both phases which are hard to uniquely dis- 
entangle. Evolution in the population mean alone does 
not necessarily imply that any individual galaxy must 
grow in size. The key evidence for growth in the "second 
phase" is the declining abundance of compact systems. 
Observationally, we seek to explain the minimum rate at 
which high-z compact galaxies must evolve so as to avoid 
leaving too many compact remnants at later times. 

To test for growth in this second phase requires the dis- 
tribution of sizes at two redshifts and the relative abun- 
dances of the progenitors and candidate descendants. We 
focus on the redshift interval z = 1 — 2 to illustrate the 



Can minor merging account for the size growth of quiescent galaxies? 



15 



method. Figure [T2] shows the cumulative distributions of 
the mass- normalized radius 7 at z = 2 and 2 = 1 using 
the fits presented in Section 3. These have been scaled 
to total number densities using the fit in Figure [11)3. We 
term these compactness functions (CFs) in analogy to 
the more familiar stellar mass function (see Bezanson et 
al. 2011 for a demonstration in terms of inferred velocity 
dispersion) . 

Mergers will shift the z — 2 CF in two ways in Fig- 
ure |12| First, galaxies will expand according to Equa- 
tion 15] which will shift the distribution rightward to- 
ward larger 7. We again assume that this growth is 
uniform and neglect stochasticity in mergers. A second, 
less important, effect is that some of these mergers will 
be among galaxies within the sample. This will reduce 
the number density of the sample over time, moving the 
size distribution parallel to the logn axis. In Section 4, 
we called these "intrasample mergers" and measured the 
fraction /is of physical secondaries they represent. The 
rate of intrasample mergers is then /is /pair /tb- Bearing 
in mind that /is/pair ~ 0.03 — 0.06 is small (Table 3]), 
we can approximate the resulting reduction in number 
density over an interval Ai by 



Alogn sa log(l - /is/pa 



(6) 



We note that the secondaries removed from the sample 
in intrasample mergers will preferentially be low mass, 
but since 7 is defined to be statistically independent of 
, it is a good approximation to shift the CF uniformly. 

The evolution of the 2 = 2 CF will thus proceed right- 
ward and slightly downward in Figure |12[ as indicated 
by the thick arrow that connects the z = 2 CF (solid 
red line) to that of its descendants (dashed). The length 
of the arrow indicates the magnitude of the size evolu- 
tion and depends on /pair a-nd t^. A plausible evolu- 
tionary path must shift the z = 2 CF to lie below that 
observed at z = 1, otherwise too many compact descen- 
dants would remain at z = 1. We use this to define a 
minimum growth rate for z = 2 compact galaxies consis- 
tent with the observed depletion in the number density 
of similarly compact systems. 

Before embarking on this task, it is necessary to define 
a minimum percentile Pmin of the z — 2 CF that we wish 
to fit within the observed z = 1 distribution. For exam- 
ple, if we require only that the largest 30% (Fmin = 0.7) 
of the z — 2 descendants fit within the z = 1 distri- 
bution, then no size growth is necessary, as Figure |12| 
shows. At the other extreme, if we require that all z — 2 
descendants are accommodated (Pmin — 0), the mini- 
mum necessary growth is approximately the same as the 
difference in the population means at the two redshifts, 
which we considered in Section 5.3. 

In practice, some intermediate Pmin must be chosen. 
Although smaller values of Pmin provide stronger con- 
straints, this must be balanced against our wish not to 
extrapolate fits of the observed 7 distribution down to ar- 
bitrarily small 7, where they are poorly constrained by 
the finite number of galaxies. In the following we conser- 
vatively set Pmin = 0.2, which is large enough that the 
empirical CFs are fairly well const rain ed at z < 2 (see 
Figure [5]). The thick arrow in Figure 12 has the minimum 
length necessary to shift the 20th percentile of the z = 2 
CF beneath the z = 1 distribution. The corresponding 



growth A log 7min can be taken as the minimum amount 
of growth necessary to sufficiently deplete the abundance 
of c ompact systems. Following the discussion in Section 
|5.3[ this minimum growth is less than the difference in 
the means of the two CFs indicated at the bottom of 
the panel. Errors on Alog7min are estimated by repeat- 
ing this calculation using many samples from the Markov 
chains used to fit the 7 distribution (Section 3). Number 
densities are also randomly perturbed fro m t he mean fit 
as illustrated by the gray band in Figure pTjp. 

Figure |12^ demonstrates that the minimum growth 
over z = 1 — 2 is Alog7min = 0.16 ± 0.03, assuming 
a size growth efficiency of a = 1.6. Throughout we 
take (/i*) = 0.39 and also set /is — 0.18 appropriate 
to z ~' 2 (see Table l3|, although the results are ex- 
tremely insensitive to this value. This minimum growth 
can now be compared to that expected from mergers 
via Equation [51 Alog7mcrg = 0.08 ± 0.02, assuming a 
short timescale of Tg = 1 Gyr. Therefore, it appears 
that only Alog7morg/Alog7min ~ 50% of the required 
growth over this interval can be attributed to mergers. 
For longer merger timescales or lower grow th efficiencies, 
this fraction would be less. In Figure [T2b we show how 
the fraction depends on Tg and a. The dotted line out- 
lines the region of likely parameters discussed in Sections 
5.1 and 5.2. We conclude that mergers alone are unlikely 
to achieve the minimum rate of expansion required be- 
tween z = 2 and z = 1, even for favorable assumptions 
regarding these theoretical parameters. 

The exercise can readily be repeated over other red- 
shift intervals. Figure [13] shows the minimum growth 
rate that progenitors at z = 2 and z = 1 must undergo 
to avoid leaving too many late compact remnants. This 
minimum rate is compared to the mean e volu tion of the 
quiescent population, introduced in Figure pTj and to the 
merger model predictions. During the period z = 1 — 2 
over which the number density is rapidly increasing, the 
population mean (black symbols) evolves more quickly 
than the minimum rate (red line). Both, however, ex- 
ceed the expected growth rate from mergers (gray band), 
even for a favorable choice of Te = 1 Gyr and a = 1.6. 
From z — 0.4 — 1 the minimum growth rate (blue) is 
only slightly less than the rate at which the population 
mean evolves, owing to the more gradual number density 
increase over this period. However, this slight decrease 
brings the required size evolution closer to the merger 
model. We conclude that mergers are roughly consistent 
with producing the more modest size evolution at z < 1, 
assuming the same favorable choices of a and Tg. 

Broadly speaking, our more elaborate model reaches a 
similar conclusion to that we inferred from a nai've con- 
sideration of the mean sizes at various redshifts (Section 
5.3, Figure 10). However, even this refined model in- 
volves some questionable assumptions. First, we have 
neglected the contribution that measurement errors make 
to the width of the observed distribution, on the grounds 
that they are expected to make a small contribution. In 
Appendix B we discuss how our results would be im- 
pacted if the true measurement errors increase rapidly 
with redshift. Second, we have assumed that the de- 
scendants of quiescent galaxies are also quiescent, but 
some systems may be r ejuvenated by seco ndary episodes 



of star formation (e.g., Treu et al. 20051. Since our re 



suits are driven by the abundance of the most compact 
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Figure 12. Left: Cumulative compactness functions at z = 2 and z = 1, representing the comoving number density of quiescent galaxies 
more compact than a given 7 based on fits described in Section 3. The thick arrow indicates the minimum size expansion necessary to 
ensure that descendants of the z = 2 population can be accommodated within the z = \ distribution. As discussed in the text, we do not 
use the lowest 20% of the 2 = 2 CF as a constraint. Bands indicate Icr confidence regions. Right: The fraction A log7nicrg/A log 7njin of 
the minimum required growth rate over 2 = 1 — 2 that is producible by mergers under various assumptions for the timescale and growth 
efficiency, with the 10 — 90% confidence intervals in parentheses. The dotted box outlines the range of likely parameters discussed in 
Sections 5.1 and 5.2. 
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Figure 13. Minimum required rates of growth for quiescent pop- 
ulations in place at 2 = 2 and 2 = 1, indicated by the dark 
gray bands (colored bands in online version), are compared to the 
growth rates of a simple merger model (light gray) discus sed in the 
text. Black diamonds reproduce mean sizes from Figure [TT^ . The 
thickness of the merger trajectories refiects the uncertainty in /jvf . 
All indicated uncertainties are \a. 

systems, which are overwhelmingly quiescent, we expect 
this to be a small effect. For example, only ~ 15% of 
z > 1 galaxies that are more compact than the median 
quiescent galaxy at the same redshift are classified as 
star-forming. Third, the lower- 2; CF in our comparisons 
applies to a constant mass threshold of logM* > 10.7. 
Since we expect the population in place at high- 2; to be 
continually growing in mass, an evolving mass threshold 
would be more appropriate. For the specific assembly 
rate IvU/M^, k, 0.03 dex/rg expected from our pair anal- 
ysis, this translates to reductions in number density of 
w 10 — 20% over the redshift intervals we considered. 
Since the last two effects are modest and oppose one an- 
other, neglecting them is justified. 
Finally, it is important to consider the stochasticity of 



the merger progress. Obviously, every galaxy cannot un- 
dergo exactly 0.16 mergers with mass ratios /x* > 0.1 per 
timescale. In reality, since the expected number of merg- 
ers per timescale is significantly less than unity, many 
galaxies will experience no such mergers over an interval 
of several Gyr. This retards the move me nt of the com- 
pact end of the distribution in Figure |12^ , leaving even 
more late compact remnants. Accounting for stochastic- 
ity would therefore only strengthen our conclusion that 
additional mechanisms are necessary to explain the rate 
of size evolution at z > 1 . 

In summary, our models for size growth via minor 
merging can reasonably account only for that observed at 
2 < 1. The faster growth rate at higher redshift remains 
difficult to explain via merging alone, even when one ac- 
counts for the rapid buildup of the quiescent population 
over the same period. 

6. DISCUSSION AND CONCLUSIONS 

Using high-quality near-infrared imaging from 
WFC3/IR taken as part of the CANDELS survey, in 
conjunction with other multi-wavclcngth data in the 
UKIRT Ultra Deep Field and GOODS-South fields, we 
have compiled a uniform sample of 935 galaxies with 
stellar masses greater than IQ^^-'' Mq and photometric 
redshifts 0.4 < z < 2.5. Within this sample, the most 
compact objects at a given redshift are those with 
quiescent stellar populations. For this subsample, the 
mean half-light radius measured at fixed stellar mass 
grows by a factor of 3.5 over this interval. The growth 
rate per unit time is noticeably quicker at early cosmic 
epochs, corresponding to the redshift range z « 1.3 — 2.5. 

We have explored the physical origin of this size growth 
in 404 quiescent galaxies over 0.4 < z < 2 by searching 
for close pairs whose photometric redshifts imply a likely 
association with their hosts. The depth of the imaging 
allows us to probe secondary companions whose stellar 
masses are only 10% of their primary hosts. Our main 
conclusion is that the delivery of stellar mass in mergers, 
estimated via the incidence of close pairs, cannot account 



Can minor merging account for the size growth of quiescent galaxies? 



17 



for more than roughly half of the minimum size growth 
that z — 2 quiescent galaxies must incur to avoid leav- 
ing a greater number of late compact remnants than is 
observed. At 2; < 1, on the other hand, mergers may 
account for most or all of the size growth rate, but only 
if a short merger timescale (~ 1 Gyr) and fairly robust 
growth efficiency {a ~ 1.6) are valid. These conclusions 
hold if the evolution of the mass-size relation is driven in 
part by the emergence of new, systematically larger qui- 
escent galaxies. If this is not the case, then the merger 
rate will fall further short of that needed to drive the 
observed size growth. 

Given the variety of theoretical and observational in- 
gredients in this analysis, it is worthwhile to review 
the assumptions underlying this conclusion. Foremost 
is the uncertainty in the merger timescale and growth 
efficiency. Most of the results in Section 5 assume opti- 
mistic values for the theoretical parameters (rg = 1 Gyr, 
a = 1.6). Furthermore, all mergers in our models are dry 
and thus provide the maximum amount of size growth, 
whereas many minor mergers at high-redshift may in fact 
involve gas-rich secondaries (Section 4.4). We also note 
that our correction for unbound projected pairs (Cmg, 
Section 5.1) is not specifically calibrated for red galax- 
ies, which are more strongly clustered, and may therefore 
understate this correction and thus overstate the merger 
rate. Altogether, it is therefore easy to argue that merg- 
ers produce less size growth than we have presented, but 
it is hard to see how the effect of mergers could be much 
larger. 

The tension at high redshift can be viewed as a conse- 
quence of the observation that the rate of size growth per 
unit time is considerably larger beyond z ~ 1.3, whereas 
the pair fraction remains nearly constant. One conceiv- 
able explanation is that the merger timescale declines 
with increasing redshift. However, current theoretical 
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Incompleteness due to 

photometric redshitt errors is a concern as higher red- 
shifts are probed, but our bes t estimates of the catas- 
trophic error rate (Section 4.11 are not high enough to 
significantly alter our conclusions. We note also that al- 
though the energy arguments discussed in Section |5.2| 
are generally applicable, the details of our framework for 
analyzing size growth are premised on spheroid-spheroid 
mergers. This is true for most other observational studies 
to date, since the theoretical framework for such mergers 
has been most extensively developed. Further studies 
of simulated spheroid-disk minor mergers, particularly 
with progenitors consistent with 2: ~ 2 observations, are 
needed to better assess the growth efficiency when the 
incoming stellar material is more loosely bound. Still, 
our pair fraction measurements imply that only about 
50 [Te/1 Gyr]~^ percent of z ~ 2 quiescent galaxies ex- 
perience any > 0.1 mergers over z = 1 — 2. This 
is likely to pose a challenge regardless of the particular 
merger physics. 

An equally important assumption is that the observed 
half-light radii are valid proxies for half-mass radii. The 
former are measured observationally, but the latter are 
relevant when considering the mass-structural changes 
caused by mergers. Although a detailed study of color 
gradients and their evolution in our CANDELS sample 
is beyond the scope of this paper, these data do confirm 



earlier studies that quiescent galaxies at z '-^ 2 typically 
display negative color gradients (i.e., are bluer on the 
outside), and that these tend to flatten at lower redshift 
( van Dokkum et al.||2010 Guo et al.|2011 Cassata et al. 
2011[ ). I'he color gradient s probably arise irom a complex 
combination of age, dust, and metallicity gradients, but 
in any case the stellar mass-to-light ratio is lower on the 
outside, so that these galaxies are more compact in mass 
than in light. If anything, we therefore expect to have 
underestimated the rate of structural change. 

Much of the early skepticism regarding the rapid size 
evolution of early-type galaxies focused on the possibil- 
ity of severe observational errors in measuring the key 
parameters of size and mass. Stellar masses could be 
overestimated by imperfect population synthesis mod- 
els, or effectiv e radii could be under estimated in shallow 
imaging (e.g., Mancini et al. 2010). Subsequent obser 



vations have weakened these claims. Although substan- 
tial unc ertainties reniain in stellar population synthesis 
models ( Muzzin et al.|2009 ), dynamical masses measured 
from absorption spec tra m moderate samples at z ^ 1.3 
( Newman et al.|2010 ) and for a fe w individual g alaxies or 
stacked spectra afT^^ 1.6 — 1.8 ( Cappellari et al. 2009} 



van de Sande et al. 2011 ) have not indicated large sys- 
tematic discrepancies with photometrically-determined 
stellar masses. Regarding size measurements, the CAN- 
DELS survey represents a major advance as it provides 
the first large space-based sample taking advantage of 
the improved depth and sampling of WFC3 relative to 
NICMOS. The radial surface brightness profile of a typ- 
ical z ^ 2 quiescent galaxy in our sample can be traced 
to ~ IRe- 

Several theorists have compared the rat e of galaxy size 
evolutio n in simulations to observations. [Hopkins et al.| 
(]2010b|, based on a suite of cosmological, hydrodynami- 
cal simulations, also conclude that mergers alone do not 
generate the entire rate of growth observed for quies- 
cent galaxies. To explain the remainder, they propose 
a combination of several physical and observational ef- 
fects. First, they assume that stellar masses and effective 
radii are over- and underestimated, respectively. Sec- 
ond, they suggest that the presence of blue cores implies 
that half-mass radii are larger than the measured half- 
light radii; as discussed above, the op posite appears more 
likely. Third, Hopkins et al.| ( 2010b) model adiabatic ex- 
pansion due to mass loss from stellar winds, but this 
ef fect alone expand s galaxies by only ~ 20%. 



Oser et al. (2011 ), on the other hand, present hydrody- 
namical "zoom" simulations in which galaxy size evolu- 
tion at z < 2 agrees well with their compilation of obser- 
vations and attribute the size expansion primarily to mi- 
nor mergers. As these authors note, one concern is that 
the absence of supernova feedback in this set of simula- 
tions enhances the stellar mass formed in low-mass halos. 
This could overstate the effectiveness of minor mergers 
by substantially increasing the stellar mass they deliver. 
As simulations and observations at z ~ 2 improve, it 
may be possible to test such effects through additional 
comparisons, such as the stellar mass-halo mass relation 
or the evolution of the stellar mass function. 

C. Nipoti et al (2011, submitted) construct a ACDM- 
based analytic framework, supported by suites of TV-body 
spheroid merger simulations, to predict the evolution 
of early- type galaxies undergoing dry mergers. Using 
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a compilation of observations of early-type galaxies at 
z — 1 — 2.5 (including this work), they conclude that 
mergers alone are not consistent with the observed rate 
of str uctural evolution a t z > 1.3. Following on earlier 
work ( Nipoti et al.| 2009"), they also find that mergers in- 
troduce too much scatter in the scaling relations at lower 
redshift unless the progenitors are finely tuned to occupy 
a very tight region in the mass-radius plane. Such fine 
tuning is not consistent with the near constancy of the 
scatter that we observe in this plane. 

Future work can extend this study in many ways. 
Imaging of the remaining CANDELS fields will allow 
possible trends of sizes and pair fractions with mass, 
redshift, and environment to be discerned more clearly, 
which may shed light on the responsible physical mecha- 
nisms. Multiplexed near-infrared spectrographs soon to 
be commissioned on 8 — 10 m telescopes will provide red- 
shifts and confirmation of the quiescent nature for larger 
samples at high redshift than has previously been possi- 
ble. This will provide an invaluable test of the photomet- 
ric redshift and star formation rate estimates on which 
the present study depends, although with current tele- 
scopes we are likely to continue to rely on photometric 
estimates for many of the faint companions. It should 
also be possible to significantly enlarge the library of dy- 
namical mass estimates, of which only a handful are cur- 
rently available for quiescent galaxies at z > 1.5, and 
thereby test the accuracy and precision of stellar mass 
estimates at higher redshifts. Spectroscopic indicators 
of maturity and recent star formation activity (Balmer 
lines, 4000 A break) may allow tests of the "minimum 
growth" hypothesis considered in this work, i.e., that 
early quiescent galaxies remain the most compact sys- 
tems in place at later epochs. If this is not the case, the 
challenge of accounting physically for the rapid growth 
of quiescent galaxies will be further heightened. 
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oiog y at z= 1 a|„g y at z= 1 a|„g ^ at z= 1 

Figure 14. Effects of measurement errors in stellar masses and radii on our conclusions. Contours show the factor by which A log 7min j the 
minimum necessary growth for quiescent systems over z = 1 — 2, would change under different estimates of the measurement errors criog-, 
at 2 = 1 and 2. Top and right axes show the corresponding errors 0"iogM» in stellar mass, assuming 10% errors in radii. Shaded regions 
indicate the corresponding confidence level (CL) at which merging alone as a viable growth mechanism is rejected. The star indicates the 
estimated errors based on the rms formal stellar mass uncertainties from SED fitting. Panels display results for several effective merger 
timescales Te and a = 1.6. As discussed in Appendix B, if the uncertainties on stellar mass (modulo IMF choice) are the same (dotted 
line) at 2 = 1 and 2 = 2, then our results are unchanged, and merging alone as a growth mechanism is rejected at > 95% CL. Only if the 
errors are much larger at 2 = 2 and the merger timescale is short is the CL reduced. 



APPENDIX 



A. THE SKEW NORMAL DISTRIBUTION 

In Section 3 we fit the distribution of sizes of quiescent galaxies to skew normal distributions that evolve with 
redshift in order to assess changes in the mean and dispersion. The skew normal distribution has the probability 
density function 

P{x) = — e / e^^dt, (Al) 



characterized by the parameters {ip,u}, s). Throughput this paper, we use a parameterization in terms of the mean 
and standard deviation a, which relate to (V', w, s) through the relations x = -ip + u}Sy/2/Tr and — a;^(l — 26"^ /n 
where 6 = s/^/l + s^. The shape parameter s relates to the skewness, and s — {) recovers a Gaussian distribution. 



B. MEASUREMENT ERRORS IN STELLAR MASSES AND RADII 

As discussed in Sections 3 and 5, the distributions in 7 that we fit and compare to merger models are dominated by 
the intrinsic variation in 7, but also include some component of scatter arising from measurement errors in the radii 
and stellar masses. Random errors in stellar mass estimates are small with good photometry (typically ~ 0.1 dex in 
this work; see also Auger et al. 2009), but systematic errors are not well understood. Comparison with independent 
dynamical mass measurements can place upper limits on the true scatter in stellar mass estimates. At z ~ 0, this 
limits the scatter to gio e-M, 5, 0-15 dex, base d on the SDSS sample described in Section 2.6, while at z ~ 1.3 the 
sample of spheroids from Newman et al. (2010| indicates a similar scatter of criogM, ^ 0.1 — 0.2 dex when spectroscopic 
redshifts are used. Analogous comparisons are currently not possible at higher redshift. For this study, as we describe 
below, the absolute uncertainties are not as important as how they may evolve with redshift. 

The main sources of systematic uncertainty include the unknown IMF and the complexities of stellar population 
synthesis models. The former is less critical for our analysis, since in our "minimum growth" test (Section 5.4) we 
are tracking the same sample of massive, quiescent galaxies, so the IMF should not cha nge. The latter uncertainty is 
likely more important, since younger populations may be systematically different (e.g., Maraston] 2005[ Conroy et al. 
|2009). A simple estimate of this effect can be obtained by comparing stellar mass estimates from the B(JU3 and CB07 
models, which differ in their treatment of the TP-AGB stars that may dominate the NIR light at ages of ~ 1 Gyr. As 
discussed in Section 3, these models predict stellar masses systematically offset by — 0.05z for our quiescent sample. 
However, the scatter between the two models actually declines slightly from ~ 0.1 dex at z = 1 to ~ 0.05 dex at z = 2. 

The main concern for our "minimum growth rate" study is the reverse: that the scatter in stellar mass measurements 
for quiescent galaxies at z = 2 is much larger than at z = 1. In this scenario, the true abundance of very compact 
z = 2 galaxies would be smaller than our fits indicate, since some are simply scattered to small 7 through random 
errors in and logM*. If the z = 1 measurement errors are comparable, then this effect approximately cancels in 
our comparison of compactness functions at the two redshifts. But if the z = 1 measurement errors are smaller than 
those at z = 2, we expect that the minimum necessary size growth over z = 1 — 2 would lessen. 

We can test the effects of redshift-dependent errors in our comparison of two compactness functions CFi and 
CF2 with estimated measurement errors criog7,i > <y\o^-y,2 by convolving CF2 by a Gaussian with dispersion a — 
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\l '^to"-! 1 ^ '''log 7 2- Having thus matched the measurement errors, we then derive the minimum necessary growth 

Alog7inin following Section 5.4. The contours in Figure [14] show how this minimum growth would change for evolution 
over z — 1 — 2 assuming various measurement errors at z = 1 and 2. The three panels consider a range of merger 
timescales spanning the range discussed in Section 5.1 and a growth efficiency of a = 1.6. As anticipated, when the 
measurement errors are nearly equal (dotted line), the derived minimum growth is not affected. When CTiog-y is much 
greater at z = 2 than at z = 1, however, the true minimum growth rate may be smaller and thus more comparable to 
the rate attainable through me rgin g. 

The shaded regions in Figure [14] display the confidence levels at which merging alone as a driver of size evolution is 
rejected. These panels indicate that the claim that merging alone is insufficient at z > 1 is seriously weakened only if 
the measurement errors are significantly larger at z = 2 than at z = 1 and the effective timescale is very short (~ 1 
Gyr). Given the other assumptions entering this exercise that are favorable for mergers (namely, that all of the most 
compact systems at z = 1 are descended from z = 2 quiescent galaxies, that all mergers are dry, and that there is no 
stochasticity in the incidence of mergers), we believe that our main results are robust. 



